


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1975 


The effects of simple boundary irregularities 
on acoustical mode propagation in ducts with 
pressure release boundaries. 


Nunes, Luis Fernando Pereira da Silva 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/20839 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist sha Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ia) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


THE EFFECTS OF SIMPLE BOUNDARY 
IRREGULARITIES ON ACOUSTICAL 
MODE PROPAGATION IN DUCTS 
WITH PRESSURE RELEASE 
BOUNDARIES 


Luis Fernando Pereira da Silva Nunes 


























. THE EFFECTS OF SIMPLE BOUNDARY 
A(T RREGULARITIES ON ACOUSTICAL MODE PROPAGATION 
IN DUCTS WITH PRESSURE RELEASE BOUNDARIES | | 
by 


Luis Fernando Pereira da Silva Nunes 


December 1975 





'Thesis Advisor: Pee COD 


pens | 











Approved for publie release; distribution unlimited. 


T17080¢ 





UNCLASSIFIED 


SECURITY CLASSIFICATION OF THIS PAGE (When Date Entered) 


REPORT DOCUMENTATION PAGE 


. REPORT NUMBER 2. GOVT ACCESSION NO. 


&. TITLE (and Subtitie) 


Bae BLftects of Simple Boundary Irregulari- 
ties on Acoustical Mode Propagation in 
Ducts with Pressure Release Boundaries 







READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


3. RECIPIENT'S CATALOG NUMBER 


5. TYPE OF REPORT & PERIOO COVERED 
Master's Thesis; 


December 1975 


6. PERFORMING ORG. REPORT NUMBER 






















8. CONTRACT OR GRANT NUMBER(Ss) 






7. AUTHOR(a) 
Luis Fernando Pereira da Silva Nunes 












10. PROGRAM ELEMENT, PROJECT, TASK 
AREA & WORK UNIT NUMBERS 







9. PERFORMING ORGANIZATION NAME AND AOORESS 


Naval Postgraduate School 
Monterey, California 93940 

















112. REPORT DATE 


December 1975 


13. NUMBER OF PAGES 


47 


[ 1S. SECURITY CLASS. (of this report) 


11. CONTROLLING OFFICE MAME AND ADDRESS 
Naval Postgraduate School 
Monterey, California 93940 

















& MONITORING AGENCY NAME & ADDRESS (If different from Controiiing Office) 





Unelassiivted 


Sa. OECL ASSIFICATION/ OCWNGRACING 
SCHEDULE 


16. DISTRIBUTION STATEMENT (of thie Report) 


‘Approved for public release; distribution unlimited. 


17. OISTRIGUTION STATEMENT (of the sbetreci entered fn Block 20, if different from Report) 


18. SUPPLEMENTARY NOTES 


118. KEY WORDS (Continue on revoree eide ff neceeeary end identify by biock number) 


Acoustical Mode Propagation 


20. ABSTRACT (Continue on revezze cide if neceseary and idantify by biock number) 
Effects of boundary irregularities on acoustical mode 


propagation in ducts, with pressure release walls are 
theoretically studied. Several particular wall irregularities 


are explored: Sinusoidal, 6-function and spatially decaying 


ee eee eee 








Bon? ass 
EDITION OF 1 HOV 6315 a eit 
DD . jan 73 1473 IS OBSOLETE UNCLASSIFIED 
(Page 1) S/N 0102-014- 6601 | ee ae 
SECURITY CLASSIFICATION OF THIS PAGE (Filion Data interest) 


Ab 





wie. 


1TY wairmyu es: 
a 


SECURITY CLASSIFICATION OF TRIS PAGE(When Deta Entered. 


=&20 


wf L J hd oD Lil iz SCHOOL 
NTEREY, CALIFORNIA 93940 
wee ee 


fo, ABSTRACT Continued) 


wall perturbations. Strong effects (resonances) appear in 


the case of the sinusoidal wall. In the other two cases, the 


resonances are not present, but there are still traveling 


disturbances whose effects can be important at great distances 


from the wall irregularity. 


mi For 473 
Bei oeg ie UNCLASSIFIED 
S/N 0102-014-6601 <= 


SECURITY CLASSIFICATION OF THIS PAGESWhenr Date Entered) 








The Effects of Simple Boundary 
Bake clad mit TeSmOnEeAGoucmIcal Mode Propagation 
in Ducts with Pressure Release Boundaries 
by 


Luis Fernando Pereira da Silva Nunes 
First Lieutenant, Portuguese “Navy 


Subomateea: In partial’ fulfaliment of the 
requirements for the degree of 


MASTER OF SCIENCE IN ENGINEERING ACOUSTICS 


from the 


NAVAL POSTGRADUATE SCHOOL 
December 1975 





NAVAL POSTGRADUATY ceuc 
MONTEREY, CALIFORNIA Fanaa 


ABSTRACT | 


Effects of boundary irregularities on acoustical mode 
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ieee DUCTION 


The theory of acoustic wave propagation in bounded 
media (waveguides and waveducts) has been treated with a 
certain intensity in the last twenty-five years. It is 
possible to classify the problems arising in these studies 
in two categories. One category involves the effects on 
acoustic propagation of small inhomogeneities and the other, 
the effects on the acoustic field of irregularities on the 
boundaries of the waveguide or waveduct. 

These two factors are the most important concerning the 
acoustic wave propagation and M.A. Isakovich [4] in 1956 
said that "such problems, besides their theoretical interest 
Gee also important from a practical point of view; thus, for 
example, in super-long-range waveduct propagation of sound 
in a laminarly~-inhomogeneous medium the effect of even small 
inhomogeneities or corrugation of the boundaries can produce 
a very great effect upon the range of propagation". 

M.A. Isakovich studied propagation in a waveduct with small 
inhomogeneities and in a waveduct with one corrugated wall, 
using the method of small perturbations. This was a natural 
extension of the problem of scattering of sound from one 
Single rough wall. 

Pema |S Moy, 1n 1956, studied the problem of 
propagation of acoustic and electromagnetic waves in ducts 


with two corrugated rigid walls. He found, in the case of 





harmonic roughness, "that waves are not scattered upstream 

if the frequency of the transmitted signal exceeds certain 
critical values related to the natural modes of the waveguide 
and the wavelengths of the wall roughness." 

C.S. Clay [6], in 1963 presented a paper where the 
effect of a slightly irregular boundary on waveguide 
propagation was studied. The results were checked by an 
experiment in shallow water and he found that the effects of 
small roughness on the surface are to increase the attenua- 
tion as a function of range and decrease the coherence of 
acoustical propagation. 

In more recent years (1971), R.H. Ebert [7] calculated 
the fluctuation of sound due to influence of a standing 
gravity wave in a waveguide. The results showed a resonance 
if the wavelength of the surface wave was one-half of the 
acoustic wavelength. In an experiment realized in the Naval 
Postgraduate School he verified the existence of that 
resonance. 

R.F. Salant [8], in 1972, stated that the research and 
study of the theory of sound propagation in waveguides could 
be and had been applied to several practical problems and 
gave as examples: transmission of noise through ventilating 
ducts, jet-engine inlets and other ducting systems. He 
considers propagation in a waveguide having two sinusoidal 
walls having a difference of phase, as opposed by earlier 
work by J.C. Samuels [5], whose waveguide had the sinusoidal 


walls in phase. 
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Salant concluded that the phase difference between the 
walls, the wall wave number and the mode number in the wave- 
guide affect strongly the disturbance generated by the walls. 

He found that the disturbances when odd modes are propa- 
gating in the waveguide are extremely different from the 
ones formed when even modes are propagating. 

A. H. Nayfek [9] in 1974 treated again the acoustic 
propagation in ducts with sinusoidal walls using the method 
of multiple scales, in contrast with previous works where 
the small perturbations method was used. For travelling 
waves, resonance occurs wherever the wall wave number is 
equal to the difference of the wave number of any two duct 
acoustic modes. 

N. E. Davis and M. W. Kenyon [1], in December 1974 
studied the effects of travelling sinusoidal surface waves 
in a waveguide, verifying that resonances occur under certain 
conditions relating wavelengths of surface and acoustic waves. 
A waveguide with pressure release walls was built and the 
presence of such resonances verified. 

In this work, a simple theory is presented for solutions 
of the acoustic velocity potential in irregular ducts with 
pressure release walls using the Fourier transform of the 


functions describing the irregularity of the wall. 


qa 





ease DeSCREP TION 


fee GENERAL FORMALISM 

Consider the existence of a waveduct, with transverse 
dimension Lo: Let the walls be pressure release and smooth, 
except for the upper wall that will be irregular. This 
waveduct will be excited by a monofrequency sound source, 
and is assumed to be of infinite length in the x direction. 

For the case of perfectly smooth pressure release walls, 
a solution for the wave equation (velocity potential) is, 


taking for simplicity an amplitude of one, 
od = sin(k, Zz) cos{wt - ne 4 ar G21) 


This represents a wave travelling in the positive x direction. 


In this equation, ky and k, must obey the relations: 


= Ww, 2 Dele 
as Vee -~ ko] 


(2628) 


Assume that the upper wall of the waveduct is irregular 


and located at 


Z= 2 + © £(x) | (2.3) 


LZ 





where f(x) is the function that describes the disturbance. 
Make e€ << 1, so that the amplitude of the disturbance will 
be small compared with the dimension ve of the waveduct. 


Pecumenaece tution Lor this case 


Cee Siete OU eet he Oe (2.4) 


In the above expression, 6 1s as given in equation (2.1), 
representing the solution for the undisturbed waveduct. 
Since the walls are pressure release, the boundary 


conditions are: 


z= 0 
¢= Q at 
Z =f + € £(x) 
Zz 
C2 4.5) 
Z= Q 
®o = Q at : 
0 - 
Z = Le 


Writing ¢ as a Taylor series expansion about the point 


1 





0 


o[2, + © £(x)] = o(2,) + cfd, + £(x) oA] 


OZ L 








Z 
2 
do 0 & 
+ e*[0, + £(x) ao + F EM (x) — AI, 
OZ Z 
+ 
: J 
rat a £4 (x) @ Pn-4 
a : 7 : !, 
3=0 OZ Z 
ieee (256) 


(eee the imposition of boundary conditions, Eq. (2.5) 


then yields the following set of equations: 


0 ® 

















es, 
| ‘ pal 7 E(x) 
2. a 
2 
3 9-6 
Oe ee! on des 0 Z 
25 = ya | 2 ao mie E(x) 
g Zz. OZ 
Zz 
; 5) 
_ 2 £4 (x) é en-3 
eee ee On 
zZ j=1 J 922 





In the assumed solution, ®% and ” must satisfy the 


homogeneous, linear wave equations 


D: 
oo= (ye - +2, 6 = 0 
st 
G i (Des; 
Do eee = 
ong = (V ~~ x) ®, = 0 
Ce oc 
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wooestetueron Of (2.1) into the first equation of (2.7) gives 


o) = - k, cos (k,, Za GOS (iit, “= WS x) Foes Ge 
Q iN 
Z 22 
a ke cos(nt) cos(wt - kK. So) pve es (2.9) 
Equation (2.9) may be generalized by making it complex. 
%) = - £(x) k. cos(nmT) er (OtoK x) CZ 230) 
Z 
such that 
= ; 2 
Re{9) 7 >, : (Zee) 
Z 
The solution will then become 
TT = Re{9,} g 272 ) 
Define: 
migée = £(x) e : C2243} 


miewansertion of definition (2.13) into Eq. (2.10) produces 


stot 
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d = -e ky cos (nm) F(x). (2. 


14) 





By use of the definition of the Fourier transform of 


“Foo bi 
PE GOe oo dx 


II 


G(s) 


and its inverse (2 wl) 


-+-0o 


lg) = = [ Gey a= 


Gx 4 


mae boundary Condition, Eq. (2.14), becomes: 


So 


ST © AL Les 
gy ; = -e k, cos (n7) in J GCS) eC ds . C2516.) 
Z, 
Use as a trial solution: 
ee ie oe, Sin(Kz) 1xs 
Epo 8 ecos(nn) 2 sin(ki.) G(s) e ds 
C2.) 


Since the above equation must satisfy Eq. (2.8) [the 
homogeneous, linear wave equation], the quantity K must obey 


the relationship 


eee ee (2.18) 
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B. SINUSOIDAL UPPER WALL 
Let the upper wall of the waveduct be perturbed in the 
fori Of a CoSinusolde of amplitude Ek where € << ] 


(See Figs. 1 and 2) The upper boundary condition is 


6= 0 at Ze & (1 + € COS Yx) 


Thus f(x) is specified as 


f(x) = Le Coslyx) C259 ) 


Mgcemeing Bq, (2.19) into Eq. (2.13) 


=k 
xX 


if 
us 
T 


Le Ces (yx) e 


igen GO (2.20) 


wie FOuUmLer transform Of (2.20) 1s: 


x oo Se =) x Saeed yo) oe 
G(x) == ; [e se +e as mew a. cud 
=F Sib [§(stk,-y) + §(stk ty) ] e221) 
and the solution will be: 
2 ee : ok "Bi Ke | _ | 
2) es ee 7: J gin Ki, {12,8 (stk, y) + 6(stk,+y) 1] 


wna dat... (2.22) 


ane, 





enh 


FIGURE 1. The Unperturbed Duct 
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oe 
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.FIGURE.2. The Perturbed Dine t 
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Using the properties of the 6-function, the integral is 


readily evaluated so that 


ate 
= —e k OSs omit 
® ZZ 


WwW, 2 22 
sin [ (=) ko 4 i [wt-(k ty) x] 


e =f 


Nl be 


sin[(2)* - (k,+y)71°/? 2, 


Ww, 2 ee 
sin[(—) = (ky ee | Z 
fo ee 


: DA me. 21/2 
sin [ (=) aes ‘oe R 


1 [wt-(k,~-y) x] 
} (202s) 


Z 


Since Op = Re {9,} , the first order correction is 


Kz 
> | See Cosine = (kL t+y) x] 1 


sin K z 
+ sin K £7 cos [wt ats (ky v) x1 } ' (2224) 


2 2 


where K* = [(2)* - (k,#y)*] (2.25) 
ioe validaty and Confirmation of Solution 
The above solution is consistent with the results 
presented by N. E. Davis and M. W. Kenyon [1]. The waveguide, 


in that case is perturbed by a traveling surface wave of 


19 





amplitude A. The upper boundary condition is 6 = 0 at 


Z= Le oe f(s where 
A 
Pate) = exe GOs (ve) and oo rt 
Z 


aemaking °& = 0 in Eq. (2.12) of {l1] and discarding the 
functional dependence in the y-direction, the result agrees 
SempolLetely with Eq. (2.22). 
2.) Properties of the Solution 
DavEGeu(2.2ee, the propagation constant in the 


Z—-direction is 


KE = ((2)? - (Kysy)?y7/? CS 
Since (2) ¢ = Iie a go F 
(Clete Oise) 
K* =k, [1 - ee (2.27) 
cea) 


Lf K* is real and: 


K7 9 = ma , £m=l,2,..., (2.28) 


mie CONGLCiONS for a resonant situation exist. From 


Boe (2.2), KS - — aan ees Inserting this into 
Z 
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Bo. (2.28) the following expression results, which states 


the relation necessary for a resonance: 


a Gia A 


2 
(ko /Y) 
(Oln & 
my 2 rs (1 + 2 k,/y) 
" (k,/y)° 


(22) 


The resonance curves are shown in Figure 3 as a plot of boy 


versus (k,/y)? . This is in accordance with Eq. 


beeere 2 Of [1], and [3]. 


aa 


(266) Manic. 
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FIGURE 3. 


he Resonance Curves 
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fil. ANALY Sis 


Using spectral analysis it is possible to find the 
solution to the problem involving any special boundary 
condition related to the upper wall of the waveduct. The 
only thing that has to be known is the Fourier transform of 
the combination 


-i1k x 
£ (3) ae ; (321) 


where f(x) represents the function describing the upper wall 
of the waveduct. 
Knowing this, the problem is reduced to solve the 


integral 


co E . 
I = f Sin Kz 1xXS 


sin Ki; G(s) Se ds 1 (3.2) 


i 
2am 
Pi GerG(s) 1S the POourier transform of (3.1), 


oo ~i(k,+s)x 
G(s) = f f(x) e 


=O 


ax ; (3.3) 


In the problem treated and solved before, the upper wall was 
described by a very simple function and the evaluation of 
the integral did not pose any special difficulties. But, if 
mx) is a more complicated function, a careful analysis of 


(3.2) is necessary in order to solve it. 


AE 





Poe line SOLUTION INTEGRAL 

Lf s is generalized to a complex quantity (s is the 
variable of integration), a careful study of the behavior of 
the function under the integral sign in (3.2) in the complex 
S-plane is necessary. 


Lic tuNnceElonm tO Study is. 


sin[(2)? - s7)°/? 
Pg) laa 
sin[ (4)? Sec oa 


1xS (3.4) 


Z 


1. The Poles 
Assuming, for the moment, that G(s) does not have 


any poles, interest is focused in the ratio of sines in (3.4). 


Since [2] 


CO 


sin z=2 JI (1 - ——=~ ) r (325) 
ee. 
m=1 m 1 


f(s) can be written as 





. [(2)*-s*] 2° | 
[ (4) 2-524272 g - h a — G(s) et*S 
£(s) = 8 AES (3.6) 
(Chee _ tess") 2, 
i @ m=1 Dee 2 
m T 


or simply 
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There are no branch cuts, the poles of f(s) are simple and 


located at 


Mo,2, 1/2 (3.8) 


It is possible to arrive at the same result by simply looking 


at the denominator of (3.4). The poles will be located at 


2 MAP 
[(M)2-5217°/% 9 = som (3.9) 
Wie ee 
Manipulation reveals that Eq. (3.9) is equivalent to 


Eq. (3.8), but here it is not evident that the poles are 


Simple, and that there are no branch cuts. 


From Eqs. (2.2) it is known that 


and 


a 





Using this in Eq. (3.8), it can be shown that the location 


of the poles are at 


os Di Wi 2 on Ne 
ea = [k. + (n“-m~) J Ca} 
fea Ms 2 eee 


Looking at this expression, it can bee seen that the poles 


will be located in the real axis if 


k > Gy ene ems) (3). 15) 


and on the imaginary axis if the inequality is reversed. 

Let us choose n=l for all that follows. This means 
that only the lowest mode can propagate energy. It can be 
seen that as m increases the associated poles are found 
successively closer to the origin on the real axis until 
above some critical value, Ms the poles lie on the imaginary 
axis with distance from the origin increasing with m. The 


first pole, Sys will be located at 


Sj = + ky e (3.12) 
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If it is assumed that eS or the fourth pole is already 
en the imaginary axis. In Table I and Figure 4 a situation 
like this is shown. 
Pee ERfectes or AbSOrption 
If the effects of absorption are introduced, the 
wave equation that governs the propagation inside the duct 
is modified and becomes: 


2 rel 


3 = 
= 32) ‘i = Q (3 7013)) 


where the a's are the coefficients of absorption responsible 

for losses. The a's can be assumed as having the same 

magnitude since all the $'s have identical form, and the 

propagation constant in the x-direction will become complex. 
It is known that the ¢'s have the form 


i fwt-k,x] 
sin(k,z)e a4) 


Beamag this in Bq. (3.13), the value of ky is 


where C2715) 


AT 





f 


TABLE I 


The Poles of 


m=1 S, = 
=? S5 = 
m= 3 S3 = 
=4 5, > 
=5 fe = 
Me OG 0s p ere.< Sas 


va: 


i+ 


i 


Fue 


‘GR 


f(s) 
k 
eG 
Ik a 3(-) 21/4 
Zz 
Ik? - a(t) 2y7/? 
S 
Ik _ i ei 
Z 
8 [15 (2) * ~~ 23272 
2 x 
ti [24( 2)? - k 2,1f2 
2 x 
* i ((m@=1) (Ey? - «29777 
Z 








FIGURE 4, 


The Poles of IE ICS) 


76,8) 





1s the total absorption coefficient in the x-direction and 
will be considered very small. 

It 1s necessary to know what is the effect of the 
introduction of absorption in the location of the poles of 
f(s), especially the ones located in the real axis. If Sy 
is located at Ss, = tk, , it can be seen that the poles 
in the negative part of the real axis will be moved up and 
the poles located in the positive part will move down 


(Figure 5). We may now choose a contour in order to solve 


the integral. 


ieee tae CONTOUR INTEGRAL 


The integral to be solved is 


oS a £(s) ds (2G) 
where f(s) is the expression in (3.4). 
foe Cholce of the Contour 
The integral in (3.16) is very similar to the one 
Isakovich used in his work [4]. It reduces to 271 times 
the summation of the residues at the poles of f(s). 
The contour to consider will consist of the segment 


[-R,+R] on the real axis and the upper semicircle Cer as 


shown in Figure 5. 


Ieee cout 
60 ap! se 
I= f f(s) ds = lim f cme ehese £(5) = 9 (3-17) 
0 R+o —-R m 
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+R 





FIGURE 5, Effects of Absorption on 


the Location 
Of the Poles of ES), 


Contour. 
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Mp li ny wd mis 4) wy detiningsih(s) as: 


sin[ (4) * - yee 
H(s) = ———_——__—__ (3718) 


sin[ (2) * - s“}L, 


HGmmreicucace where sGtc)mdees net have any poles the 


integral will be: 


co ixs_ 
i i f(s) ds 2 ) G(s, ) e 
m 


—©co 


Res {H(s_) } GS ot) 


Sieene Other hand, 1f the Fourier transform of F(x) has poles 


of its own, at s = Ba! where 
Le (3.20) 
Sts : 


the expression to evaluate the integral is modified as 
follows: 
co ixs 


T= f Cs) sas = 27 2 G(s.) e a Res{H(s,) } 


1xs 
+ 27i ) H(s_) e 4 Res{G(s_)} (eal) 
g = 99 


In either case, the poles of interest will be the ones 


included within the contour chosen. 
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Poe DeTinteten Of the Transrorm of F(x) 

In order always to be able to close the contour in 
the upper half plane, the definition of the Fourier transform 
ue ix) must be altered. 

We define 
sal feo’ 


G(s) = f F(x) e dx 


a OO 
with inverse 


F(x) == f G(s) e°*** ds 


The upper sign of the exponent will represent the 
positive x-part of the waveduct with waves standing or 
propagating to the right and the lower sign will represent 


the negative x-part of the waveduct with waves standing or 


propagating to the left. 
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EV. DPE CIAL CAS ES 


In order to simplify the problem and apply this theory 
to some specific cases we shall impose some conditions on 
the problem. Since 

en , (4.1) 
in the duct, the cut-off frequency will be 


Ww = ck : (4.2) 


The sound source that is going to excite the duct will have 


a monofrequency of 


w= Ac kK, ' (4.3) 
where (Pecwieo< 2 (4.4) 
to assure propagation only in the lowest modes. Under this 


condition and knowing that 


ie eee (4.5) 


(only lowest mode allowed in z, n=1) the location of the 


poles of f(s) will be at 
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bet W, 2 Mneee 7 2 
ee t= 
= + (eee ace (4.5) 


me 
ua 
Z 


It can be seen that only one pole will be on the real 
axis (sj), and that all others are already on the imaginary 


axis. 
The general expression for the location of the 


poles of f(s), assuming G(s) has no poles is: 


oman aye am) (4.6) 


m x 
do 


and the meh residue is 


leA2 
sin{ t(2)7-s7] 2? G{s) ets 
Res_ = 


m 1/2 -1/2 
2 2 can Hy Qs 2 2 
cos} (2) SS ] Le } a {Ug S ] Le } 28h, 





(4. 7} 
evaluated at s = Sa 
This will become, after manipulation, 
mk ie 
ee oh eee G(s.) e-em (4.8) 
m Q F —*-m 
m Z Z. 
Meee lec ys 
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Mie eirece son acoustic velocity potential, due to the at? 


pole .es [using Eq.o(2.17),)J 
Lwt 


e k > Pay Sie , (4.9) 


a 
Zz 27 be 
Pie cusimlcas) Poles Ol eits Own. (4-2)ewillsbe modified 


and the effect on the velocity acoustic potential will be 


Lwt ae 
e k, a7 2 1a p Res, 


+1ixs 
+ Y H(s,) e J res,{ (G(s, 3} ], 
(4.10) 


Peele o-FUNCTION DISTURBANCE 


Define f(x), the function describing the wall disturbance 


as 
F(x) = =. On s2) : (4.11) 
Then 
“ik, x 
F(x) = =e 6(x) e (4.12) 
and 
Cle “2, (4513) 


with no poles. 
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The poles of f(s) will be located according to Eq. (4.6). 


The first pole, the only one in the real axis will be at 


= -—k (4.14) 


and the second one, on the imaginary axis will be at 


s, = +i My 2" -1) - k7)°/? (4.15) 
= + iP 
where P is defined as 
2 = (so? = 2 Aes (4.16) 
Le Xx 


The overall situation is shown in Figures 6 and 7. The 
residue of f(s) at Sy: Res, (s,) Waetlesoen USang Ho. (4. 0)5 
kK, Fik,x 
Res, (sj) = — sin k 2 eG (ary) 


k 
X 


The acoustic potential due to this residue will be calculated 


Mem (4.9), and 


k i (wt: k_x) 
b.4R et =: 4 
cy es, (sj) es ae Seen k 2 ie : (4.18) 
X 


oy 





i 


ELCGURH@G eee O-Function Disturbance 





{ 


FIGURE 7. The Poles of Interest 
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Since o, = Re{o,} , in this case 


¥ 2 
— _ Z . . = 
o, = a sin k 2 Sin (wttk,x) (4.19) 


The above equation represents two scattered traveling waves 
emanating from x=0, one going in the positive x-direction and 
the other in the negative x-direction. 


The residue of f(s) at Sor Res, (So), will be, using 





fore (4.8), 
ae ce ocie 
Res, (so) = -i —5- sin 2k 2 e (4.20) 
with a corresponding acoustic velocity potential of 
2 ; SPN Go ee) 
@, [Res (s9)] = sin 2k,z e” e* (4.21) 


This represents two disturbances in the acoustic field, 

decaying in the +x-direction. It can be seen that far away 

from the wall disturbance their effect will be negligible. 
Very close to the wall disturbance, x=0, the effects will 

be considerable and depend on the size of Ee diminishing 

as P., increases with m. Ee 1s defined as 


2 = ye 2272 


X t 4 M= oo 5. s 


= 2 T 
a [(m i 


However, even at great distances from the irregularities 


in the wall the acoustic field is affected by >, [Res, (s,) 1. 
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The total field, discarding the effects of decaying 


disturbances as in (4.21) is: 


rk Z 


Sin kzlcos (wt-k,x) - € 3 sin (wt7k,x) ] (4.22) 


mecestiat the traveling perturbation is 90° out of phase with 


the initial wave. 


B. THE DECAYING WALL 


Define f(x), the function describing the wall disturbance 


as 

f(x) = -2, e” [ax | (4) 
Then 

ex! pce 

F(x) = -2, e |*I ¢ (4.24) 

and 
2h, 
Ge) = 6 === (4.25) 


Jaftz+ (2+ 14] 


with poles located at 


s = a(-k, + 1) (4.26) 
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The poles of interest of f(s) will be located at 


Ss, = =i 
So = +iP (;.27) 
Ss = a(-k, +1) 

The overall situation is shown in Figures 8 and 9. Since 


ere ehas poles, Bq. (3.21) must be used to evaluate the 
integral. 


The residue of f(s) at Sis Res, (s,) will be 


Res, (s,) = 2 2 fle) a ae as (4.28) 
Mie eam oe oD 2 a 
a + IS (a-1) 


aw 
No 


* 


with a corresponding acoustic velocity potential of: 


k Z 
©, [Res i.) or 2 
I t= Ky a? -|- k* (a1) 


BANC rel OS), 
: x 
sin k,2 e 


(A227) 


representing two traveling waves emanating from x=0, going 


in opposite directions, and 


2k * 
= 





o — Re{9, } =- 7 sin kK, Z Sin(wtyk, Xx) 


va) 
(4.30) 
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FIGURE 8. The Decaying Wall 





FIGURE 9. The Poles of Interest 
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mne residue Of £ts)> at Sor Res, (s5) will be: 


GL 


: 2 
iPa[l a (a + ky) 


sin 2k 2 aus (4223) 





Res, (s,) s 


with a corresponding acoustic velocity potential of 


2 
4k, Tit Fee 
@, [Res,(s,)} = —————+,——_———— sin 2k_z e e 
1 eee Pafl + (>=) f Kk) 4] Z 
(4.32) 


This represents two standing disturbances decaying with 
x and negligible at great distances from the origin. 


The residue of f(s) at Ss Will be computed using (3.21) 


and 1S 
+ixs 
H(s_) e 3 res{G(s_)} Ea 2s) 
g g 
or 
Sin YB+tiC 2 *#ax 7K x 
_— nn Ok e 1 L. (4.34) 
Sin YB+iC L,, 
where: 
be 23 2 a eee 2 
Eee a (1 eae) eet! 
(4.35) 
- 2 
C = 2a k 


43 





The corresponding acoustic velocity potential is: 


Sin /B+iC gz are i (wt Fk, x) 
o, [Res, (s,) } =k ——_—_-—-———- e* “ e (4.36) 


ge Sin YB+iC Le 


This corresponds to two traveling-decaying waves (both 
meme and Z Girections) going 1n Opposite directions from 
x=0. At large distances from the origin its effect will be 
meecrically nil. 

Comolning resules, Only valid at x 2>0 or x << Og; 


we have: 
od = 29 7c aT 


ae 
sin k_z [cos(wt-k_x) - € a 
& i Ky a® +k“ (a-1)4 


a 





sin(wt7k,) ] 


It 


(4.37) 


It can be noticed, again, that the traveling disturbance 


is 90° out of phase with the original propagating wave. 
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V. RESULTS AND CONCLUSIONS 


The results obtained show that the effects of very small 
irregularities in the physical boundaries of the duct can 
have appreciable importance in the waveduct sound propagation 
even at very large distances from the disturbances. In the 
case of the sinusoidal wall, resonances can exist, which 
could have the same order of magnitude of the original 
propagating wave. These resonances do not appear in the 
two other special cases studied, the é6-function disturbance 
and the decaying wall. It seems that the theory works and 
can be applied to other types of corrugations, once the 
Fourier transformation of the function describing the wall 
disturbance is known. However, complications can arise 
depending on the complexity of the above mentioned Fourier 


meeansctOorm, making the solution integral difficult to evaluate. 
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